ac Josephson effect in the resonant tunneling through mesoscopic 


superconducting junctions 


Yu Zhu', Qing-feng Sun?, Wei Lit, Zhong-shui Mat, and Tsung-han Lin* 
1 State Key Laboratory for Mesoscopic Physics and 
Department of Physics, Peking University, Beijing 
100871, China 
2 Center for the Physics of Materials and 
Department of Physics, McGill University, Montreal, 


PQ, Canada H3A 2T8 


Abstract 


We investigate ac Josephson effect in the resonant tunneling through meso- 
scopic superconducting junctions. In the presence of microwave irradiation, 
we show that the trajectory of multiple Andreev reflections can be closed 
by emitting or absorbing photons. Consequently, photon-assisted Andreev 
states are formed and play the role of carrying supercurrent. On the Shapiro 
steps, dc component appears when the resonant level Eo is near the position 
VL + iimw = VR + imaw (e = h = 1), where mı and mz are integers, w is 
the frequency of microwave, Vz and Vr the chemical potentials of the super- 
conducting leads. Analytical result is derived in the limit A — oo, based on 


which new features of ac Josephson effect are revealed. 


PACS numbers: 74.50.+r, 73.63.Kv, 73.23.Ad 


Introduction-Since the achievement of atom-size superconducting quantum point con- 
tact (SQPC) by break junction technique [1-4], considerable theoretical works have been 
devoted to the transport problem in S-I-S or S-N-S structures (S=superconductor, I=tunnel 
barrier, N=normal metal). Comparison between experimental data and theoretical calcula- 
tion suggests that the coherent picture of multiple Andreev reflection (MAR) is the central 
conception in understanding the transport through such mesoscopic SQPC. Nevertheless, 
less attention is paid to the resonant tunneling (especially in ac transport) through meso- 
scopic superconducting junctions, in which the central region is consist of discrete electronic 
states. Such resonant tunneling can be achieved in nanoparticles sandwiched between su- 
perconducting membranes [5], or metallic clusters absorbed on the STM substrate [6], or a 
piece of carbon nanotube tunnel coupled to superconducting electrodes [7], or quantum dots 
fabricated in the S/2DEG hybrid structures [8], etc. 

The transport through these structures is greatly modified due to the existence of discrete 
energy levels in the central region. It is shown by Yeyati et al. [9] and Johansson et al. 
[10] that the combination of MAR processes and resonant transmission gives rise to a rich 
subharmonic gap structure (SGS). Let us consider the transport through a resonant level 
with width I’ coupled to S leads with superconducting gap A. In the limit T > A, the 
resonance is sufficient broadened such that the I-V characteristics mimic that of SQPC. In the 
regime I ~ A, pronounced SGS appears with concomitant negative differential conductance, 
which is distinguished from the nonresonant transport through SQPC. In the limit T < 
A, although SGS is more particular in the log-scale plot, the tunnel current as a whole 
is exponentially small. This is because a narrow resonance can not cover all the MAR 
trajectory, and it is very unlikely that multiple discrete levels are just in the position where 
MAR trajectory passes (see Fig.1a). 

In this paper, we address the effect of microwave (MW) irradiation on the resonant 
tunneling through mesoscopic superconducting junctions. We show that the tunnel current 
is greatly enhanced on the Shapiro steps, when the resonant level Ep is near the position 
Vt + iimw = VR + imaw (e = ħ = 1 throughout), where mı and mg are integers, w is the 
frequency of MW, Vz and Vp are the chemical potentials of left and right S lead. This can 
be attributed to the formation of photon-assisted Andreev states (PAAS), which play the 
role of carrying supercurrent. In the limit A — oo, analytical results are derived, revealing 
new features of ac Josephson effect in the case of resonant tunneling. 


Hamiltonian and Formulation—The model Hamiltonian reads H = X` B=L.R Hg + Heent + 


Hr, in which Hg = $ po Erab Asko + Yea ct M + H.c.) is the Hamiltonian for 
the th S lead, Heent = Eo >. CÈ Co the resonant level in the central region, and Hr = 
J sko (Ua EXP i Í Valea | Beads + H.c.) the tunnel couplings. Vg(t) = Vg + Wg coswt is 
the time-dependent voltage drop across the th barrier, where Vg is the chemical potential 
controlled by dc bias voltage and Wg coswt the ac voltage induced by MW irradiation. For 
simplicity, we assume that Coulomb interaction and multi-level effect can be ignored in the 
central region. This assumption is somehow too ideal to achieve, but allows us to obtain 
analytical results and therefore instructive for the understanding of more complicated cases. 


Keldysh Green functions are defined in the 2 x 2 Nambu representation: 
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The time-dependent current flowing out of the th lead can be expressed as 
Ig(t) = 2ReTra,(G" o X5 + Gs o X3 )i, (2) 


in which g, is the 3rd Pauli matrix, o the shorthand notation for integration over intermedi- 
ate time variable, and Ng the self energy due to tunnel coupling between the central region 


and the Gth lead. G", G*, and GS satisfy the integral equation: 
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with Ig being the coupling strength, 7 the dephasing rate in the S lead, f(e) the Fermi 
distribution function, and Im yZ > 0 as a convention. The remaining task is to solve these 

integral equations and evaluate the dc component of the time-dependent current. 
There are two intrinsic frequencies in the problem, wı = 2V = 2(V;, — Vr) and w = w. 


Generally, one may perform a Fourier transform 


Ann = Jo dtas f Tet Afe) (11) 
lily 

and derive the recursive relation for Aj,),(€), as done in [11]. The calculation in this way 
relies on the numerical computing power, and analytical results are impossible. We note, 
however, that in the case of narrow resonance, dc component appears only if w, = Nw», with 
N being an integer. This becomes clear by considering PAAS shown in Fig.1b and Ic: an 
electron (a hole) incomes through the resonant level Eo and Andreev reflected by the right 
superconducting gap, the reflected hole (the reflected electron) exchanges mı photons with 
the MW field so that it can again passes through Æo. Then the hole (the electron) is Andreev 
reflected by the right superconducting gap as an electron (a hole), and exchanges m2 photons 
to close the trajectory (detailed discussion for photon-assisted Andreev reflection is available 


in [12]). It is easy to see that the formation of PAAS requires 2V = (mı + m2)w. For this 


reason, we shall only consider the case of V = Vy = Kw, while the current deviated from 
this condition is negligible small. 
The problem is largely simplified since w can be used as the basic frequency in the Fourier 


expansion. Define the Fourier transformation as 
A ilwtı de —ie(tı—t2) A 
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The definition guarantees the nice property that if C = Ao B then Cy,(e) = 
Xr Am (€)Ben(€). The Fourier transformed g”* and 53°* are 
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in which ag = We is the MW strength on the th tunnel barrier, J;,(x) the nth Bessel 


function, ep =e-—lwand ¢ = e — lw + Vn, Vi = Vw and Vp = 0 are set as a convention. 


Correspondingly, the equations for G”, G°, and G'S are Fourier transformed into G™“(¢€) = 


g™"(e) + g(x (e) G (e) and GS (e) = G" (e) ¥<(e)G° (e), or equivalently, 
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We note that finite order perturbation expansion is inadequate in the problem, because 
the formation of PAAS involves up to infinite order of tunneling processes. To re-sum up 


the series, we adopt the resonant approximation [12] 
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The approximation implies that the overlap between sidebands Eo + lw and Eo + lw can be 


(19) 


ignored if l Æ ly, which is justified when Ts < w. Applying the approximation to Eq.(17) 


and Eq.(18), one can obtain the solution 
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It can be shown that the relations G” = (G*)' and G< = — (G<)! still hold in the solution. 
With these Green functions, the time-dependent current can be expressed in a summation 


over ac components 
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The current formula can be applied to the study of ac harmonics. However, we are more 
interested in the dc component J = I; = —Ip, due to the experimental reasons. To produce 
analytical results, we take the limit A — oo, due to which all single particle processes are 


forbidden and Andreev reflection is the only conducting mechanism. After some algebra and 


take the limit 7 — 0 (notice that lim, wiz = 76(x) ), one can derive 
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in which E$ = +4/ (Eo — 2mw)* + |[Lm|? and Lm = 4 [T rIm(2aR) + Pre Jmn (2ar)]. 
Eq. (26) is the central result of this paper, which is for the dc component of Josephson 
current in the resonant tunneling through mesoscopic superconducting junctions. 

Before numerical study, we make a few remarks on this result: (1) The phase dependence 


of I is mainly determined by the prefactor sin @, and a weak cos dependence is hidden in 


Ex, via I'm. For this reason, ¢ is set as 5 in the numerical calculation. (2) For the special 
case N = 0, ag = ag = 0, and T; = Tr =T, I = 2sin gee lf (—£0) - f (+£)], 
with Ey = Ee + (T cos T which reproduces the exact result for de Josephson current 
through a resonant level. (3) The current formula is for the gauge choice of V, = Vy and 
Vr = 0. It is easy to see that the formula is invariant under the transformation L — R, 
Vy > —Vy, N — —N, and Ey — Eo—Vy, which corresponds to the gauge choice V; = 0 and 
Vr = —Vy. (4) E= are the poles of G” (e), which can be interpreted as PAAS. Obviously, E= 


m m 


contribute to the supercurrent with opposite signs. One can expect that resonant structures 
will appear near Ey = imu. (5) Bessel functions enter the current formula not in the 
square form J2(x). It can be shown that I(az,aR,¢) = I(—az,—ar, ġ) for N even and 


Iar, QR, $) I(—az, — QR, Q + T) for N odd. 
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numerical results-Firstly, we discuss the case that MW is applied symmetrically to the 
left and right tunnel barrier, i.e., a, = ag. Fig.2 shows the curves of I, = I(¢ = 5) vs 
Eo — ix at bias voltage Vy = xu, with N = 0,1,2,3,4 (the curves are shifted by = 
for demonstration). Four features are noteworthy in the plot: (1) The curve is symmetric 
(anti-symmetric) with respect to Ey = 4" = for N even (odd). (2) Photon-assisted 
structures appear near Ey = Vr — $(N — m)w = Ve + mw. (3) There are two types of 
resonant structures, single peak (dip) and peak-dip pair. (4) The structures grow with the 
MW strength non-monotonously.. Feature (1) is due to the relation (Eo, ¢) = I(Vn — Eo, @) 
for N even and I(£o,¢) = I(Vyn — Eo,¢+7) for N odd. Using I(Eo, —¢) = —I(Eo, $), one 
can obtain I(Eo,¢ = $) = (—)*I(Vy — Ep, = 3). Feature (2) can be understood in terms 
of PAAS shown in Fig.1b and 1c: Electron and hole are Andreev reflected back and forth 
by the superconducting gaps. When Ep is near the position Vg — E(N — m)w = Vpr + iimw, 
quasiparticles may exchange N — m photons with the MW field at the left tunnel barrier 
and exchange m photons at the right, so that the trajectory is closed and bound states 
are formed. It is these PAAS that play the role of carrying supercurrent. Understanding 
of feature (3) and (4) needs quantitative analysis of each photon-assisted structure. To 
proceed, we decompose the total current into I = >>, Im, and expand Im near Ey = iimw. 
Let Ey = 3mw + 6, one can obtain Im = I} + Ip, with 

1 

4 (82 + Eml’) 
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Oi corresponding to single peak (dip) and peak-dip pair. At zero temperature, Ci 
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The cancelation between J* and J> results in two types of resonant structures: C4 


and Cə are proportional to 
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The oscillatory nature of Bessel functions in these coefficients is responsible for feature (4). 

There is an interesting special case, N = 2, m = 1, a, = ap, due to which Ci = Cy = 0 

for arbitrary MW strengths. Correspondingly, the resonant structure near Ep — Yy = 0 is 


always missing for N = 2. 


Next, we investigate the case that MW is applied only to one of the tunnel barriers, i.e., 
ay = 0 and ag # 0. Fig.3 shows the curves of Ie vs Eo for ay = 0 and ag = 1. In contrast 
to the symmetric case, a single peak is pinned at Ey = Vz. The reason is as follows: since 
no MW is applied to the left barrier, Andreev tunneling through this barrier occurs only 
when Ło lines up with Vz, while photon-assisted processes is allowed at the right barrier 
irradiated by MW field. The inset shows Ie vs MW strength arg at Ey = Vz for N > 0 
(The case of N < 0 can be easily deduced from N > 0). One can see in the plot that the 
peak height (including the sign) is predominated by the prefactor Jn_—n(2az)Jm(2aR). For 
ar = 0, JIm—n(2az) requires m = N, and the peak height is proportionate to Jy(2aR). We 
note that there exist some regions of MW strength where Ie and Vy have opposite signs, 
which can be viewed as quantum pump effect. However, this feature is dramatically different 
from the quantum pump effect in normal mesoscopic junctions. 

conclusion—To sum up, we have investigated ac Josephson effect in the resonant tunneling 
through mesoscopic superconducting junctions. We show that PAAS play an essential role 
for conducting supercurrent through a narrow resonance, in which MAR trajectory is closed 
by exchanging photons with MW field. On the Shapiro steps 2V = Nw, dc component 
appears when the resonant level Ep is near Vz — (N —m)w = Vr + mw, due to the formation 
of PAAS. In the limit A — oo, analytical result, Eq.(26), is derived for the dc component of 
Josephson current, and help to understand new features of ac Josephson effect in the case of 
resonant tunneling. We have dropped Coulomb interaction in the calculation, the derived 
results are meaningful in the following senses: (1) They are directly applicable to the systems 
where Coulomb blockade effect is negligible, i.e., U ~ I < A (U is the charging energy). 
This is possible by using substrate with large dielectric constant to reduce U or using proper 
material as S lead to obtain large A. (2) They are instructive for more complicated cases. 
We note that the conception of PAAS is also useful for the interacting case. For U > I, 
the resonant level Ep is effectively split into two resonances Ey and Eọ + U. Similar to 
Fig.1b and 1c, one can draw diagrams of closed MAR trajectory through these resonances. 
Moreover, Coulomb blockade can be partially removed by bias voltage or MW irradiation, 
as long as U is comparable to A. Obviously, calculation including interaction term is much 
more difficult, and analytical results are almost impossible. Efforts along this line are still 
in progress. 
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FIGURE CAPTIONS 


Fig. 1 Schematic diagram of the resonant tunneling through mesoscopic superconducting 
junctions. (a): Without MW irradiation, MAR trajectory can not be covered by a 
narrow resonance, and the tunnel current is exponentially small. (b) and (c): In the 
presence of MW irradiation, MAR trajectory can be closed by emitting or absorbing 


photons. Two types of PAAS are formed, carrying supercurrent with opposite signs. 


Fig. 2 The curves of I. = I(¢ = į) vs Ey — Yy at bias voltage Vy = Sw (N from 0 to 4), 
with symmetric MW strengths a; = ar. Parameters are: w = 1, Ty, = Tr = 0.02, 


kgT = 0.001. V; = Vy and Vr = 0 are set as a convention. 


10 


201708.00160v1 


chinaXiv 


Fig. 3 The curves of I. = I(¢ = Z) vs Eo at bias voltage Vy = Zw (N from —4 to 4), with 
MW strengths az = 0 and ag = 1. Other parameters are the same as Fig.2. The inset 
shows Je at Eo = Vz vs the MW strength ar. Diagrams for corresponding PAAS is 


also shown in the plot. 
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